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1 Introduction 



The affine Toda field theories associated to an afHne Kac-Moody algebra g have 
received a great deal of attention in recent years. The classical equation of motion 
which describes them is 



□ 2 0+^]>>^e^)=O, (1.1) 

where <j> is the dynamical field, taking values in the Cartan subalgebra of g (for 
whom Chevalley generators are the Hi and the Coxeter number is ft,), yu and f3 are 
real parameters, and on are the simple roots of g. Taking f3 to be real we obtain 
a model with a particle spectrum whose properties are algebraic in origin (see, for 
example, refs. 0,0]), and whose S-matrix is exactly known || ^, ^, [|, [| [?|, ||, ^, [Tc| ]. 

If the coupling constant (3 is pure imaginary then the model possesses a series of 
vacuum solutions. We expect soliton solutions to be those of minimal energy which 
interpolate these vacua, the difference between the final and the initial vacua being 
known as the topological charge of the soliton. The seminal calculation was made 
by Hollowood [|ll] for the A r models using Hirota's method. This technique was 
later developed to include the rest of the theories JlJ, |l3). A new method, involving 
the use of the general solution discovered by Leznov and Saveliev (see e.g. 0j), has 
recently proved to be very useful in extracting particular properties of the solitons 
|lj|. This method can be motivated by consideration of the form of the canonical 
energy-momentum tensor 

T, v = (d^, - 9 -f (d a <l>, 9» -^$> fc (e"°*<« - l) . (1.2) 

" k=0 

Using arguments given in ref. [jl5| it can be shown that 7) t „ can be split into a 
sum of two parts, T^ v = C M „ + Q^, where 8 M „ is traceless and C M „ is a total 
derivative. Due to the topological nature of solitons, their properties, in particular 
their energies and momenta, depend only on their behaviour at infinity, and hence 
we might expect M j, to vanish for soliton solutions. Further work |l^, [l6| shows 
that 0^,, can be written purely in terms of the chiral fields which arc the parameters 
of the Leznov-Saveliev solution, so it is plausible that the soliton solutions arise if 
these fields vanish. 

In order to find the N-soliton solutions a further ansatz for the form of a constant 



of integration appearing in the general solution must be made [17 ; the N-soliton 
specialised solution turns out to be 



(A k \ex V (w 1 F^{z 1 )y..ex V (w N F^(z N )^ \A k ). (1.3) 



1 



For a detailed explanation of the quantities appearing in the formula the reader 
should consult ref . JljJ ; we will need only the following details. The parameters ij 
are integers lying between 1 and r which describe the species of the N solitons. Zj 
are complex parameters related to the soliton rapidities pj in the following way: 

k + l 

*i=0je- pi -T; 1-4 
% 

where 9j — ±1, and the (ft. are structure constants of g which can be calculated 
from the following formula of ref. [Q 

q+ = lp • = 2ix p {l)e- 5 " Bl ^l h (1.5) 

(i = y/—T). The functions Wj encode the space-time dependence of the solution, 
and are given by the formula 

Wj = Qj exp |\/2^j | qf. \ (t sink pj — zcoshpj) j . (1.6) 

The first major use to which this formalism was put was a calculation of the 
masses of single solitons of species p. These turned out to be jl^| : 



±V2h 



An interesting fact which is not yet fully understood is that these masses are in a 
certain sense dual to those of the particles. 

In this letter we will use further features of the solution ([O]) to calculate the 
Poisson brackets on the N-soliton phase space and will perform a canonical quanti- 
sation to extract the S-matrix. 

2 Poisson Brackets 

The classical phase space of the affine Toda theories has been discussed recently in 
ref. p8[ . The symplectic form Q is calculated by integrating the symplectic current 
over all space at some time, and is simply 

fl = J dx{8(t>,5dt<t>) (2.1) 

where ( , ) denotes the usual Killing form on the algebra g. 

We wish to investigate the phase space of the soliton solutions. To do this, 
we insert the soliton solutions into the symplectic form VL. First let us calculate 
the symplectic form for a one-soliton solution. We will need to use two important 
properties of this solution, but our method will mean that we do not need the explicit 
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form of the solution itself. This is a typical feature of calculations concerning solitons 
in afhne Toda theory. The first feature we need is that the soliton is a relativistic 
object, and so the Poincare algebra must be realised on the phase space. What this 
means for the solution is that the field <f> must appear only as a function of u, where 

u = t sinh p — (x — xq) cosh p. (2-2) 

The rapidity p and centre of mass xo are the real parameters of the soliton solution. 
The relationship of these parameters to those in the algebraic ansatz for the soliton 
solution, as well as the explicit dependence of tf> upon u, can be determined from 
the formulae of jl5|, [|l7), but we shall not need these yet. Using the antisymmetry 
properties of the wedge product we obtain 



O = dx [ — , — coshpfeo A (5(sinhp), 
\du du J 

dx(8 t ^) 6X0 A5(8inhp) 



sinh p 

f 5xq A 6 {sinh p) 
~ ~J Ttx shrnT ' (2 ' 3) 

where the last step is performed using the expression for the energy-momentum 



tensor fll.2j). Finally, substituting P = — J dxT i = Msinhp into (2.3), we obtain 



SI = 5£ A Sp, (2.4) 

where £ = Mxq coshp is the canonical variable conjugate to p. This is of course the 
result we expect. 

In order to evaluate the symplectic form in the more general case of N solitons 
we consider the form of the solution as t — > ±oo. In the generic situation the 
solitons will be well separated and since 11 is a local expression we can just add up 
the contributions from each of the solitons in turn. Thus we find 

N 



We now need to relate these variables to those which parameterise the solution. 



Let us consider the form of the solution (1.3) as x increases from — oo. The 
first significant departure from the vacuum will occur with the soliton of greatest 
rapidity, p^. We know from ref. |l7| that the greatest non- vanishing power of F lN 
within a representation of level x is x.^ This in turn means that the solution for 
the component field (f>k will be the logarithm of a polynomial of degree in Wn 



defined by equation (1.6). As we move through this soliton Wn becomes much 



We are only considering the theories where g is simply-laced from now on. 
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greater than 1 and so we can ignore all but the highest power in the polynomial as 
far as calculating the form of the solution for greater x is concerned. This term of 
course multiplies an algebraic factor F lN (zn). Normal ordering the vertex operator 
expression |l7j for this yields 

e-W* = (A fc |F l "(z Ar )|A fc )^(A fe |exp(x jl , lJV (z 1 ,z Ar )iy 1 F il (^i)) ... (2.6) 
. . .exp (x N - ltN {z N - 1 ,z N )W N -iP ilf - x (z N -if) |A fe >, 

where 

h-l , x w n (7i)-73 

X tJ ( Zl ,z 2 ) = IJ 1 - e- 2 ™ 1 /^ . (2.7) 



n=0 



Zl 



The roots ji are the simple roots at up to a sign p7[ , and w is the Coxeter element 
of the Weyl group. Thus, aside from a constant shift in the field (j>, the only effect 
is to change each of the Qj for j < N by a factor Xi-i N (zj , z/v). Extending this 
argument as we move through each of the solitons in turn leads us to conclude 
that as t — * — oo the field <f> becomes a sum of one-soliton solutions but with new 
parameters 

Q? = Qsl[X j , P iz j ,z p ). (2.8) 

P>3 

The rapidities remain unchanged under the transformation to 'in' variables. Now 



all we need to do is relate the variables Qj to the £j. Comparing expressions (1.6) 
and ( |2.2[ ) we find that up to an irrelevant constant 

\/2fi\q^\ coshp 



Using the mass formula equation (1.7) we obtain 



and so 



, 2h9hxQ 

Z= , (2-10) 



# = Pi- (2.H) 



Similar arguments yield the following relationships for the out variables 



pr - Pr (2-12) 
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Since the symplectic form f2 is trivial to invert in terms of terms of either 'in' 
or 'out' variables the appropriate Poisson brackets for the natural variables follow 
straightforwardly, and can be seen to be consistent if we note that X(zj, z p ) depends 
only on the rapidity difference pj — p p . 

3 Scattering Matrix 

Having found canonical coordinates on the classical phase space we can now quantise 
the theory in a straightforward manner, simply by replacing the canonical Poisson 
bracket with the canonical commutators 

[r,p in ] = [c ou v ut ] =%n. (3.1) 

We can then attempt to discover the unitary transformation (analogue of the canon- 
ical transformation) which will be the S-matrix for the reduced theory. From equa- 



tions ( 2.11 ) and ( 2.12 ) we conclude that 



P<j \ P>3 

Pj = Pj = Pj- (3-2) 

We can see from these equations that the S matrix S is purely a function of the 
rapidity differences, and satisfies the equation 

dCLogS) —2hi I . i , \ -» . \ 

q = Tgvlfr 1^ lllX pd(Pp> Pj) ~ 2^ Pp) ■ ( 3 - 3 ) 

P j IPI \ P <j P >j J 

The solution of this is 

TTmfr £ £ w "(7a) • (lb) Li 2 exp (p a - Pb + j- (S ib , B - Si a ,B - 2 

b>an=0 J 

(3.4) 

We have used the definition of the classical Euler dilogarithm 



oo s 



^i^dy. (3.5) 

V 



Many of the interesting properties of this function were studied by one William 
Spence in the early nineteenth century [ p0| . A more recent discussion in the math- 
ematical literature can be found in rcf. pM. Recent developments relating to 



conformal field theory are reviewed in ref. (£2|. The matrix (3.4) satisfies various 
conditions, which are essentially related to properties of the time-delay functions 
Xij. This function was studied recently in ref. [^3), where it was noted that it 
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has properties reminiscent of those of S matrices. Indeed, one can think of these 
properties as being consequences of the fact that the S matrix (3.4) is periodic, 
symmetric, etc. Let us define 



h-i 

r„6(p) = 5>"(7«) • (7b) Lh 

n=0 

Then the matrix function T ab {p) satisfies the following relations: 

(i) T ab (p + 2m) = T ab (p) 

(ii) T ab (p) = T ba (p) 

(iii) T ab (p + in) = -T ab (p) 

(iv) (T ab (p*))* =T ab (p) 

(v) T ab (p) = -T ab (-p) + 2T ab {0) 

(vi) J2t=i j k Tit(p + iTjt) — 0, where I is a free label and k satisfy a fusing rule; 

■qt — — 2£t + c ^ 2 ~ 1 , where c(t) is 1(— 1) if the root a t is black(white), and £t 
are the integers in the fusing relation Ylt=% j k w_ ^*7t = 0- 

These relations can be shown directly: Property (i) is obvious; properties (ii), (iii) 
and (iv) follow using arguments analogous to those used for the time-delay functions 
X ab (p) in rcf. p3| ]. Property (vi) similarly follows from an argument analogous to 
that given in ref. when discussing the affine Toda particle S-matrix. Property 
(v) is most easily proved by noting firstly that the p derivative of this equation is 
true - this follows using eqn. ( |3.5| ) and the fact that X ab (p) = X ab (—p) (see ref. 
p!|). Thus the left-hand side of (v) equals the right-hand side up to an additive 
constant, and putting p — one sees that this constant must vanish. A direct proof 
of (v) seems more involved - for example, for the A\ case this property reduces to 
the equation 

Li 2 {-e p ) - Li 2 {e p ) = -Li 2 {- e - p ) + Li 2 ( e - p ) + 2Li 2 (-l) - 2ii 2 (l). (3.7) 
The Euler dilogarithm satisfies the following 'inversion' relation |2^] 

Li 2 (-y) + Li 2 (-l/y) =~~ \^vf- (3-8) 

The function Li 2 (y) is divergent for real y, y > 1; however, one can define a con- 
tinuous function on the real line by setting ]2^ | 

Li 2 (y) = y - Li 2 (l/y) - i(logy) 2 , for y > 1. (3.9) 



exp p 



it: f c(a) — c(b) 
~h 



2n 



(3.6) 



G 



Then, using the above two equations and the facts that Lz2(l) = 7r 2 /6 and 1) 
— 7r 2 /12, equation (3.7) may be proved. 



The constant term in the relation (v) means that the S matrix (3.4) will be 
invariant under p — > — p only if it is normalised so that S(p = 0) = 1. Physically 
this requirement is obvious - that there be no scattering when the two solitons 
have the same rapidity. We note that the properties (i)-(vi) above reflect identities 
satisfied by the Euler dilogarithm. 



Thus the expression (3.4) satisfies relations expected for an S'-matrix. Note, 
however, that the Euler dilogarithm can be continued to a multi- valued function on 
the complex plane, minus the segment (l,oo) of the real axis pl| , p2[ . Hence our 
proposed S'-matrix does not have the expected pole structure. We will comment 
upon this in the following section. 

4 Conclusions and Developments 

It is important to be clear about what this S-matrix is, and what it is not. We have 
only made a semi-classical approximation to the quantum theory of affine Toda 
solitons, and so do not expect to see all of the behaviour typical of a quantum field 
theory. In particular as noted above, there are no poles in formula (|3.4|) , and we 
have made no mention of the renormalisation of the soliton masses. The interested 
reader should consult the recent papers and |2^] for a discussion of this latter 
point. Very little is known about the expected behaviour of the poles on account 
of difficulties concerning the unitarity of the affine Toda theories in the imaginary 
coupling regime. Another feature expected of the quantum theory which this kind 
of approximation will not possess is that there will be no processes which change the 
topological charges of the scattering solitons. This is because such these are absent 
classically. An obvious conjecture worthy of exploration is that a full quantum 



version of (3.4) involves the quantum dilogarithm of ref. |2q| . 

What we believe we have is the S-matrix of an quantum integrable particle 
theory, generalising that of Ruijsenaars [^7], ^|. In the case of A-series solitons 
of equal mass this particle model already reproduces the appropriate scattering 
shifts. In the more general case however, no such model is known. Starting from 
the formula for the shifts which follows from Ruijsenaars model, and fixing the 
rapidities in such a way that the soliton fusing rule ]l7"t is satisfied, we have obtained 
correct expressions for the shifts when two solitons of different masses scatter (it is 
reasonably easy to show that all of the A-series solitons can be obtained by repeated 
fusing from the lightest in this manner). The question is whether this procedure 
can be implemented at the level of the Ruijsenaars Hamiltonian. This is a non- 
trivial problem since the fusings correspond to imposing imaginary constraints on 
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the particle rapidities, and it is unclear what effect this will have on the symplectic 
structure of the phase space. We remark that the change p — > p + 2imr/h is a 
symplectic transformation, albeit a complex one, which gives us some hope that the 
above procedure can be made to work. 

A natural generalisation of this work would seem to be to try and include the 
breathers. In the sine-Gordon theory at least these breathers are related to the 
particles, and possess a number of discrete energy levels according to the value of 
the coupling constant j3. It is a long-standing problem in affine Toda soliton theory 
to try and characterise the breather spectrum for more general theories, and in 
particular to see if the particle-soliton correspondence remains. We have obtained 
preliminary results in this direction, including a description of the energy levels of 
the A-series breathers of remarkable and rather mysterious simplicity. We plan to 
discuss these issues in a forthcoming paper [ p9[ . 
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